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Abstract 

Room temperature ferromagnetism has been one of the most sought after topics in today's 
emerging field of spintronics. It is strongly believed that defect- and inhomogeneity- free sample 
growth should be the optimal route for achieving room-temperature ferromagnetism and huge 
efforts are made in order to grow samples as "clean" as possible. However, until now, in the dilute 
regime it has been difficult to obtain Curie temperatures larger than that measured in well annealed 
samples of (Ga,Mn)As (~190 K for 12% doping). In the present work, we propose an innovative 
path to room temperature ferromagnetism in diluted magnetic semiconductors. We theoretically 
show that even a very small concentration of nanoscale inhomogeneities can lead to a tremendous 
boost of the critical temperatures: up to a 1600% increase compared to the homogeneous case. In 
addition to a very detailed analysis, we also give a plausible explanation for the wide variation of 
the critical temperatures observed in (Ga,Mn)N and provide a better understanding of the likely 
origin of very high Curie temperatures measured occasionally in some cases. The colossal increase 
of the ordering temperatures by nanoscale cluster inclusions should open up a new direction toward 
the synthesis of materials relevant for spintronic functionalities. 

PACS numbers: 75.50.Pp, 75.30.Kz, 75.40.-s 
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I. INTRODUCTION 



The hope of attaining ferromagnetic order at room temperature and above, has spurred 
a huge interest in the field of diluted magnetic semiconductors (DMSs)P® and diluted mag- 
netic oxides (DMOs) 4 6 . Extensive experimental as well as theoretical efforts have been made 
to predict high Curie temperatures (Tc) i n these materials. Among the various materials 
widely studied, one of particular interest is (Ga,Mn)N, a wide-band-gap DMS. Different ex- 
perimental results have reported Tc's varying as widely as 10 to 940 KP^with a typical Mn 
content between 7 and 9%. However, recent theoretical studies, based on model calculations, 
have predicted a T c of 30 K in homogeneously diluted and uncompensated Gai-^Mn^N for 
x— O.Od 111 , which is in good agreement with results obtained from ab initio based studies^ 
combined with the self-consistent local random-phase approximation (SC-LRPA) method. 
On the other hand, using the same ab initio couplings, the Monte Carlo studies lead to Tc 
of 35 K for 2=0.0d-HH5] in Gai-^Mn^N. These theoretical calculations predict the highest 
reachable Tc in homogeneously diluted Gai-^Mn^N. Then how can we explain the very high 
Curie temperatures observed by some experimental groups? We will provide an answer in 
the following. 

From these observations crucial questions arise. How do we explain the huge fluctua- 
tions of the critical temperatures in these materials? Is there a systematic way to boost 
the critical temperatures beyond that expected in the homogeneous compounds (inhomo- 
geneity free)? After observation of ferromagnetic order in Mn-doped Germanium (Tc=116 
K for 2=0.035)-^, several experimental studies reported quite high critical temperatures in 
(Ge,Mn) film^HSffl However, the underlying reasons were not really clear. In Ref.^ scanning 
photoelectron microscopy measurements revealed stripe-shaped Mn rich microstructures 
which were believed to be the origin of ferromagnetism in Ge^^M^. More recent experi- 
mental studies have revealed self-organized Mn rich nanocolumns formation in Gei- x Mn x , 
which gave rise to a very high Tc (>400 K)P^for x=0.06. Magnetotransport measurements, 
in this case, have also shown a large anomalous Hall effect up to room temperature. The 
spinodal decomposition (alternating regions of low and high concentration of magnetic im- 
purities) was suggested to be the reason for the high temperature ferromagnetism in this 
case. Similar nanometer-sized clusters, with increased Mn content compared to the sur- 
rounding matrix, were also detected by transmission electron microscopy (TEM) analysis 
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in Geo.95Mn .o^- In recent experimental studies on (Zn,Co)CpH, the authors claimed the 
existence of two types of nanosized ferromagnetic Co clusters. The first were spherical with 
diameters of about 5 nm leading to critical temperatures of ~100 K and the others were 
columnar about 4 nm wide, with a maximum height of 60 nm, leading to significantly larger 
critical temperatures of ~300 K. These results were confirmed by high-resolution transmis- 
sion electron microscopy (HRTEM). Hence this kind of anisotropic nanoscale inhomogeneity 
can lead to interesting magnetic and transport properties. Inspite of the existence of several 
experimental studies, the effect of impurity clustering on magnetism in DMSs and DMOs 
has been weakly studied on the theoretical front. Ab initio based studies for these type of in- 
homogeneous disordered systems are difficult due to the large size of supercells required and 
no standard methods have been proposed as yet. In Ref.^the authors have simulated the 
spinodal decomposition in DMS by using Monte Carlo methods and they have predicted an 
above-room-temperature T^ for the spinodal phase in (Ga,Mn)As and (Ga,Mn)N, calculated 
from the "standard" random-phase approximation (RPA). Here "standard" means that the 
crucial self consistency was not implemented in the RPA calculations. However these high 
Tcs were found for samples containing a relatively high concentration of Mn, above 20%, 
far from the dilute regime. On the other hand in the dilute case, for approximately 5% of 
Mn, the authors have found a suppression of the critical temperatures in the presence of 
spinodal decomposition phases. Note also that the calculations were limited to small system 
sizes compared to the typical size of the inhomogeneities and the average was done over few 
configurations only (typically 10). In RefP^ the authors have presented density- functional 
theory (DFT) based calculations of N-doped Mn clusters, and have given a hypothesis that 
a high Curie temperature detected in some of the GaMnN samples is a result of the forma- 
tion of small Mn clusters carrying giant magnetic moments. The large variation in Curie 
temperatures could be attributed to the formation of N induced Mn clusters of different 
sizes in samples grown under different conditions. Their analysis suggests the importance of 
the growth mechanism in these kind of materials. Similar density-functional calculations on 
the effect of microscopic Mn clustering on the Curie temperatures of (Ga,Mn)N were also 
reported in Ref.^. However the T^'s were calculated from the mean-field approximation, 
which is already known to overestimate the Curie temperatures in homogeneously diluted 
semiconductors. 

In this article, we present a generalized and comprehensive study of the effect of nanoscale 
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FIG. 1: (Color online) Snapshots corresponding to four different concentrations of nanospheres x n s 
(a) 0.02, (b) 0.04, (c) 0.06 and (d) 0.08. The grey (red) atoms denote the impurities inside (outside) 
the nanospheres. Here L=36, vo=2a (a is the lattice spacing) and Xi n =0.8. 

inhomogeneities on the Curie temperatures in diluted magnetic systems. The calculations are 
performed on very large systems (finite size analysis is provided) and a systematic sampling 
is done over several hundreds of disorder configurations. In contrast to previous studies, we 
report giant effects on the Tc in dilute materials. In some particular cases the Tc can be 
enhanced by up to 1600% compared to that of the homogeneously diluted system. There are 
several factors that lead to these effects, such as the concentration of inhomogeneities in the 
system, the size of the inhomogeneities, the concentration of magnetic impurities inside the 
inhomogeneities and also the range of the exchange interactions between the impurities. In 
the following we shall see how these physical parameters play an important role and affect 
the critical temperatures. 
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II. MODEL AND METHOD 



For simplicity we have assumed here a simple cubic crystalline structure and the con- 
clusions that will be drawn will be general. The sizes vary from L=32 to L=44. The 
inhomogeneities considered here are of spherical shape of radii ro- For the sake of clarity 
and to avoid additional parameters, in our calculations the total concentration of impuri- 
ties in the whole system is fixed to x=0.07. In the following we denote the concentration 
of nanospheres by x ns =N,s/N, where N5 is the total number of sites included in all the 
nanospheres and N=L 3 is the total number of sites. The concentration of impurities inside 
each nanosphere is defined by Xj n . We denote the total number of impurities and the number 
of impurities inside the nanospheres by N*^ p and N-^ p (=Xi n Ns) respectively. We choose 
the nanospheres in such a manner so as to restrict their overlap with each other. 

In Fig. [T] four typical random configurations corresponding to four different concentra- 
tion of nanospheres x ns (0.02, 0.04, 0.06 and 0.08) are depicted. As we increase the number 
of nanospheres in the system x ns increases and consequently the concentration outside de- 
creases, since the total concentration (x) is fixed. Now to evaluate the Tq, the effective 
diluted Heisenberg Hamiltonian H.Heis=-J2i j Jjj S; ■ Sj, is treated within the SC-LRPA the- 
ory. The self-consistent local RPA is a semi-analytical approach based on finite temperature 
Green's functions. It is essentially an extension of the standard RPA to the case of disor- 
dered systems. Here the thermal fluctuations are treated within the RPA and the disorder is 
treated exactly without any approximations. The Curie temperature of a system containing 
N imp localized spins is obtained from the expression 



We define the retarded Green's function as G i j{tjj)=^ oo Gij{t)e lult dt={{Sf]Sj )). The 
set of parameters \=\mi T ^ Tc (S-)/m, where m is the average magnetization, are calculated 



configuration and then averaged over a few hundred configurations of disorder. The accuracy 
and reliability of the SC-LRPA to treat disorder and/or dilution has been demonstrated 




(1) 



where 




(2) 



self-consistently (more details can be found in Refl^). The Tc is calculated for each random 
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TABLE I: The ratio, (R hom ), of the homo geneous Curie temperatures for different x to that of 
£=0.07, for X=a and a/2 



x (R hom ) (A=o) (R hom ) (A=o/2) 
0.8 9.7 22 
0.7 8.9 20 
0.6 7.9 17 
0.4 5.2 11 



several times in the past) 1 * 12 * 2 -^. 

The exchange couplings in a DMS, as found from ab initio based calculations, are rel- 
atively short range in nature and almost exponentially decaying } 14 ' 15 1 . Thus in the present 
study we have assumed generalized couplings of the form Jjj=J exp(— | r | /A), where 
r=ivrj and A is the damping parameter. In (Ga,Mn)As, for about 5% Mn a fit of the ab 
initio magnetic couplings provides a value of A of the order of a/2. Note that in the case of 
(Ga,Mn)N the ab initio couplings are of even shorter range. Thus we focus here on two par- 
ticular cases, A=a and A=a/2, where a is the lattice spacing. Although these length scales 
are comparable, in the presence of inhomogeneities the effects on the critical temperatures 
will be very drastic. In order to measure directly the effects of nanoscale inhomogeneities, 
the averaged Curie temperatures (T^ n/l ) are scaled with respect to the averaged Curie tem- 
peratures of the homogeneously diluted system (T^ om ) for x=0.07, their ratio is denoted by 
(Rc)- The averaged Curie temperatures (T^ om ) for the homogeneous systems are found to 
be 0.9 J and 0.05 J , for X=a and a/2 respectively, for x=0.07. 

In Table ([!]) we provide the averaged Curie temperatures for 80%, 70%, 60% and 40% 
homogeneously distributed magnetic impurities scaled with respect to that of the 7% ho- 
mogeneous case, for X=a and a/2. The ratio is denoted by (K hom ). These values will be 
relevant in the discussions to follow, where we consider these types of concentrations inside 
the nanospheres. 
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for X=a. (a) Results 



for a fixed concentration inside the nanospheres (xj n =0.8) and different radii ro- The upper x 
axis represents the values of x ns corresponding to Xi n =0.8. The long black arrow indicates the 
direction of increasing ro. (b) Results for a fixed radius (ro=2a) and different concentration inside 
the nanospheres. The long black arrow indicates the direction of increasing Xi n . The solid red line 
indicates the T^ CA scaled with respect to (T^; om ). In the figures, squares correspond to L=32 and 
circles to L=36. 

III. RESULTS AND DISCUSSION 



Figure 



2 a 



shows (Rc) as a function of N*^ p /N*^ p corresponding to the case of A=a. The 
concentration inside the nanospheres is fixed at Xj n =0.8 and the Tc is calculated for spheres 
of different radii. For this concentration inside Xi n , each nanosphere contains 5, 15, 26, 45, 
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and 64 impurities for r =a, v2a, 2a, \/5a, and a/6o, respectively. N^ p /N'^ p =0 corresponds 
to the homogeneously diluted case (absence of inhomogeneities). We observe a clear increase 
in the critical temperatures with increasing fraction of impurities inside the nanospheres 
as well as with the nanospheres' size. For about 80% of the total impurities inside the 
nanospheres, Tc is enhanced by up to 150% for the smallest nanospheres with r =a, and 
by almost 350% for the ones of radius io=^Qa, which is rather significant. This increase for 
io=\/6a is more than one-third of that found for the 80% homogeneously distributed case 
(Table [T]). Thus the clustering of magnetic impurities does lead to a considerable increase 
of the critical temperatures due to the strong interactions within the nanospheres. The 
other important point to take note of is the Tc obtained from the mean field virtual crystal 
approximation (VCA), Tq A —^x J^njJj, where n^ is the number of atoms in the i-th shell. It 
is well known that the VCA overestimates the true critical temperatures, often very strongly. 
However, the present results show that in the presence of inhomogeneities the VCA value 
can no longer serve as an upper bound. Indeed, as can be seen here, for a relatively small 
concentration of nanospheres (x ns ~0.2) the VCA value is already exceeded and for higher 
density of nanospheres the VCA actually strongly underestimates the critical temperatures 
in these systems. 

Let us now focus on the case where the nanospheres are of fixed radius (ro=2a) and 
the concentration inside the nanospheres vary (Fig. [2(b)] ) .(R c ) is plotted function of 
for different X{ n . The curves show an overall monotonous increase with increasing 
concentration of nanospheres. However, the enhancement of the critical temperatures is 
also controlled by the concentration of impurities inside the nanospheres. Decreasing the 
concentration inside the nanospheres effectively means reducing the number of impurities 
inside a cluster of the same size, and thus reducing inter nanosphere interactions. This could 
explain the relatively small increase in the Tc values with decreasing xi n . However, as will 
be seen in the following, the variation of the critical temperatures is more complex than 
this simple picture. Thus we find that not only the relative number of impurities inside the 
nanospheres but also the concentration inside the nanospheres have a drastic effect on the 
critical temperatures in these systems. 

Now we move to the case of the shorter-ranged couplings, A=a/2, which will appear even 



more interesting and which lead to unexpected effects. Figure 3(a) shows the (Rc) as a 
function of N-^ p /N*^ p for a fixed x in =0.8. T c is calculated for nanospheres of different radii 
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FIG. 3: (Color online) (R c ) = (T | om) as a function of N*™ p /N^ p =^x ns corresponding to X=a/2. 
(a) Results for a fixed concentration inside the nanospheres (xj„=0.8) and different radii ro- The 
upper x axis represents the values of x ns corresponding to Xj n =0.8. The long black arrow indicates 
the direction of increasing ro- (b)Results for a fixed radius (ro=2a) and different concentration 
inside the nanospheres. The long black arrow indicates the direction of increasing Xi n . The solid 
red line indicates the Tq CA scaled with respect to (T^ om ). In the figures, squares correspond to 
L=32, circles to L=36, triangles to L=40, and diamonds to L=44. 

(r =a, y/2a, 2a, yha and y/6a). We have considered system sizes varying from L=32 to L=44 
to check for the finite-size effects. The L=44 systems typically contain ~6000 impurities. 
In contrast to the case of \=a discussed above, the variation of Tc with /N*^ is not 
monotonous anymore. Here we see a colossal effect of the size of the nanospheres on the 
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Tc- For the smallest nanospheres (io=a) there is hardly any noticeable effect, with the 
critical temperatures remaining close to that of the homogeneous case. Now as we increase 
the radius of the nanospheres for a given concentration of nanospheres, there is a sharp and 
strong increase in the Tc values. As can be seen, even for a reasonably small concentration 
of nanospheres (x ns ~0.2) we obtain a remarkable jump of almost 900% for iQ=2a and even 
1600% for ro=v6a, compared to that of the homogeneous case. This gigantic increase in 
the presence of nanospheres with r =v / 6& is more than 70% when compared to the Tc 
of the 80% homogeneous case (Table [T]), which is rather extraordinary. This implies that 
in materials such as (Ga,Mn)N, where the exchange interactions are really short ranged, it 
would be possible to reach Tc>500 K (as Tc for homogeneously diluted Gai-^Mn^N is 40 K 
for x=0.07^HHI) by inducting nanoscale inhomogeneities. The presence of such nanoclusters 
may also explain the very high Tcs observed in Gai-^Mn^N by some experimental groups^. 
It should be of great interest to analyse experimentally the effect of such nanoclusters on 
the critical temperatures in these kind of materials. Here again the mean field VCA is found 
to strongly underestimate the Tc for most cases. This is expected since the mean field 
VCA treatment is unable to capture all the relevant physical effects in both homogeneously 
disordered as well as inhomogeneous systems. Thus it becomes clear that in systems with 
relatively short- ranged couplings the size of the inhomogeneities plays a very important role 
in controlling the critical temperatures. The nonmonotonous behavior observed here implies 
that several physical parameters are in competition (length scales and relevant couplings). 
Thus, we cannot explain this variation by assuming the inhomogeneities to behave as "super- 
spins" only. 



In Fig. 3(b) we consider the case of nanospheres of fixed radii ro=2a, which is particularly 
interesting. (Rc) is shown as a function of N^ p /N'^ p for different Xi n (0.8, 0.7, 0.6 and 0.4). 
For a fixed Xi n , we observe a gradual increase in the critical temperatures with increasing 
concentration of nanospheres, and then it decreases as x ns increases further. In contrast to 
the case of X=a, there is a clear maximum in the Tc around N*™ p /N*^ p ~0.2 for Xi n =80% and 
70%. For this value of N|JL,/N*^L, as we increase the concentration inside the nanospheres 
we observe a huge jump in the critical temperatures, from a small increase for 2j n =40% 
to almost 900% for x in =80%, compared to that of the homogeneous case. It should be 
noted that for /N*^ =0.9 and x in =80% the increase is reduced to about 600%, which 
is still considerably large. However for £ iri =40% we hardly obtain any significant increase 
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compared to that of the homogeneous case. Hence in this case the concentration inside the 
nanospheres is found to have a crucial effect on determining the critical temperatures of the 
system. A careful statistical analysis reveals that the case of r =2a for A=a/2 is particularly 
intriguing. As will be seen, the analysis of the Tc distributions exhibits interesting features. 
In the following we provide a more detailed study for this particular case and try to analyse 
the reasons for the origin of this kind of behavior. 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 
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FIG. 4: (Color online)Normalized Rc* distributions for four different x ns : (a ) 0.02, (b) 0.04, (c) 
0.06, and (d) 0.08 corresponding to X=a/2. Here ro=2a, rcj n =0.8 and L=32. The red dashed lines 
indicate the (Re) values, which we have shown in Fig |3(b)[ The solid black lines are a guide to 
the eye. 

In Fig. [4] we show the normalized distributions corresponding to the case of r =2a, 
Xj„=0.8 and X=a/2. The distributions are obtained using a sampling over a few hundred 
configurations of disorder (~600). As can be seen from the figure, we obtain very interesting 
wide distributions for the different concentrations of nanospheres (x ns =0.02, 0.04, 0.06 and 
0.08). For x ns =0.02 we observe a kind of bimodal distribution: one peak with a large weight 
at high T c (Tg 9/l ~ll(T£ om )) values and another one at lower T c (T^~2(T^ om )) with a 
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much smaller weight. When increasing x ns to 0.04, the width of the distribution is almost 
unaffected, but we notice a clear transfer of weight from the high Tc values to the lower one. 
By further increasing x ns to 0.06, the transfer of weight increases further: the low-Tc region 
has a significantly higher weight. Finally for relatively high x ns (~0.08), the weight is now 
concentrated around the lower Tc values and the distribution exhibits a tail-like structure at 
higher critical temperatures. This transfer of weight is the reason for the maximum in the Tc 



observed in Fig. 3(a) The origin of this kind of distribution is not very clear at first. However 
the analysis of the configurations reveals an interesting feature. We have considered two 
different kind of configurations. The first set of configurations of nanospheres corresponds to 
the situation where the distance between the nanospheres is restricted to small separations. 
The second kind corresponds to large separations between the nanospheres. First, it is found 
that in both cases the distribution of Tc is relatively narrow and unimodal. However, in the 
first case the Tc distribution is centered around T l g w , whereas in the second case it is centered 
around T c i9/ \ It is surprising and counter-intuitive that the largest Tc's are obtained from 
the configurations where the inter nanosphere couplings are weaker. This is a clear indication 
that several length scales are competing. Now the nature of distributions shown in Fig. [4] 
can be explained as follows. In the case of low concentration of nanospheres (Fig. [4^ a)) 
the probability of finding the nanospheres relatively far apart from each other is relatively 
high, and conversely the probability of finding them close to each other is relatively small. 
Thus this leads, in the distribution, to a significant weight around the high Tc values. As 
we gradually increase x ns , the probability of finding configurations with the nanospheres at 
relatively large separation decreases, while the probability corresponding to small separation 
increases. As a consequence the weight in the distribution around T c ifl/l decreases and that 
corresponding to the low Tc increases, as observed in Figs. (4fb) and (4Tc). Finally for the 
largest x ns (~0.08) the weight is mainly concentrated around the low-Tc region (Fig. |4^d)). 
Interestingly, this kind of behavior is not observed in the case of X—a, for nanospheres of 
radii varying from iQ=a to \/6a. For this case (longer range) the distribution of the critical 
temperatures is always narrow and unimodal, thus all the configurations (nanospheres far 
apart or close to each other) lead to similar values of the critical temperature. This confirms 
the idea that several length scales and typical couplings compete to give rise to this rich and 
new physics. 

Let us now discuss some experimental consequences. We have shown that in systems 
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with effective short-ranged exchange interactions it is possible to obtain two different critical 
temperatures depending on the size and concentration of the inhomogeneities and also on 
the typical separation between them. For example, in the case of A=a/2 for nanospheres of 
radii ro=2a, x ns =0.02 and Xi n =0.S, the T c * 9/l value is almost five times that of the T 1 q" value. 
This could explain the wide range of Tc values observed experimentally for materials such as 
(Ga,Mn)]>pHn] anc i the apparent dissension between theoretical predictions and experimental 
observations for these kinds of materials. In this context, it should be noted that Li et alW^ 
proposed two different ordering temperatures in Gei.^Mn^., T c and T* c with Tc^CT c . The 
higher critical temperature T c is associated with the ferromagnetic ordering temperature 
within isolated spin clusters and the onset of global ferromagnetism only occurs at Tc. For 
£=0.05 the values of Tc and T c were found to be 12 and 112 K, respectively. However, 
detailed experimental studies in this direction would help to fully confirm this picture. 



IV. CONCLUSION 



In conclusion, we have presented a detailed study of the effect of nanoscale inhomo- 
geneities on the critical temperatures in diluted magnetic systems. We have shown that 
for materials with effective short ranged exchange interactions it is indeed possible to go 
beyond room temperature ferromagnetism by inducting nanoscale clusters of magnetic im- 
purities. A gigantic increase in the critical temperatures of up to 1600%, compared to that 
of the homogeneously diluted case, is obtained in certain cases. We also provide a plausible 
explanation for the wide variation of the Tc's, observed experimentally, in some materials 
such as (Ga,Mn)N. A meticulous study revealed that the relative separation between the 
inhomogeneities can play a decisive role in controlling the Curie temperatures. In some 
cases uniform distribution of nanospheres is found to favor very high critical temperatures. 
This fact could be further corroborated by detailed experimental studies. If, by controlling 
the growth conditions, the formation of the nanoscale inhomogeneities can be manipulated, 
it will open up the possibility of studying these disordered inhomogeneous systems in more 
detail. We believe that our study will pave the way for a better understanding of the origin 
and control of high-temperature ferromagnetism in dilute magnetic systems, which can serve 
as building blocks for potential future spintronic devices. 
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